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Abstract 

This  paper  explains  some  details  related  with  the  choice  of  the  surface  or  mortar  basis  functions 
for  the  Trefftz  finite  elements.  It  is  shown  that  appropriate  normalization  of  the  surface  basis 
functions  lead  to  unambiguous  definition  of  the  generalized  admittance  matrix  descriptor  of  the 
element  and  to  a  possibility  to  define  scattering  matrix  descriptor  of  the  element.  A  brick  element 
is  investigated  to  illustrate  the  approach. 

Introduction 

The  Trefftz  finite  elements  have  been  recently  introduced  into  computational  electromagnetics 
[  1  ]-[5] .  The  procedure  of  the  element  building  is  relatively  simple.  Solutions  of  the 
homogeneous  Maxwell’s  equations  are  used  as  the  basis  functions  to  expand  the  interior  element 
fields  into  Fourier  series.  Finite  number  of  plane  waves  propagating  in  different  directions  may 
be  chosen  as  such  basis  functions.  To  build  a  descriptor  of  the  element,  an  additional  set  of  basis 
functions  has  to  be  defined  on  the  surface  of  the  element.  The  purpose  of  the  surface  basis 
functions  is  to  control  energy  transmission  through  the  surface  of  the  element.  Projections  of  the 
interior  fields  on  the  surface  basis  functions  link  the  interior  expansion  with  the  exterior  basis 
functions  and  finally  lead  to  admittance  or  scattering  matrix  descriptors  of  the  element.  To 
control  the  energy  transmission  through  the  element  surface,  the  surface  functions  have  to  satisfy 
some  orthogonality  conditions.  This  was  not  covered  in  the  papers  [ l]-[5]  and  can  be  ambiguous 
even  in  the  case  of  a  brick  for  example.  Only  Trefftz  elements  built  with  Galerkin’s  projectors 
are  considered  in  this  paper  because  of  they  are  conservative  by  construction. 


Trefftz  element  surface  field  expansion 

Let’s  consider  a  polytope  element  Cl  e  M3  with  surface  S  and  Nf  number  of  faces.  Power 

transferred  through  the  surface  of  the  element  can  be  expressed  through  components  of  the 
electric  and  magnetic  fields  tangential  to  the  element  surface: 

P  =  j>(ESlxHlt)-nsds  (1) 

where  ns  is  a  normal  vector  directed  inside  the  element,  and  ESt  and  HSl  are  components  of  the 


element  interior  field  defined  on  the  surface  of  the  element  and  tangential  to  the  surface.  ESt  and 
HSt  are  expanded  into  Fourier  series  inside  the  element  according  to  procedure  described  in  [  1  ]- 
[5].  To  expand  element  surface  field  into  a  Fourier  series,  a  finite  set  of  vector  functions  is  used: 


•f'l(m) 


Kl(m) 
K^Km)  J 


l  =  m  =  \,Ml 


(2) 


where  M,  is  the  number  of  basis  functions  at  the  face  number  / ,  n,  is  the  normal  vector  to  the 


face  number  l  directed  inside  the  element,  el(m)  is  electric  field  component  and  hl{m)  is  magnetic 
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field  component  of  the  basis  function.  Minimal  number  of  the  basis  functions  per  face  for  a  3- 
dimensional  problem  is  two,  to  account  two  orthogonal  polarizations  of  electric  field  on  the  face. 
For  a  problem  with  just  one  polarization  of  electric  field,  one  basis  function  per  face  is  the 
minimal  number.  To  construct  conservative  Trefftz  finite  elements  (TFE),  it  is  necessary  to 
control  energy  transmission  through  the  surface  of  the  element.  In  other  words,  the  vector  basis 
functions  (2)  have  to  be  an  orthogonal  set  with  a  power  scalar  product  defined  on  the  surface  of 
the  element  as  follows: 

(.//(m)’  fk(n))  =  X  \(n)  )  '  nsds  =  d l(m),k(n )  (^) 

S 

where  Sl(m)  k(n)  is  Kronecker  symbol  equal  either  to  the  unit  if  l  =  k  and  m  =  n  or  to  zero 

otherwise.  If  condition  (3)  is  satisfied,  the  electric  and  magnetic  fields  on  the  surface  of  the 
element  can  be  expanded  into  Fourier  series: 

Nf  M, 


ESI  —  y,  y,  Vl(m\ei(m\ 

1= 1  m= 1 

Nf  M, 

(4) 

HSt  =  y, 

l  =  \  jYi— 1 

The  coefficients  of  expansion  v/(m 

)  and  il(m)  can  be  found  as  follows: 

F  ( m ) 

f  X  Km)  )  ■  Thds  =  f  Est  ■  Km)ds 

SI  SI 

1*  /  jjj  - -  \  1*  - - 

(5) 

h(m)  =  _ 

)  [%n)  xHS<)-  Wlds  =  J  H  St-  hl(m)ds 

SI  SI 


With  the  nonnalization  (3),  power  transferred  through  the  surface  of  the  element  (1)  can  be 
expressed  as  follows: 

Nf  M, 

^=ZEw*(«)  (6) 

l  m 

This  definition  of  power  corresponds  to  the  well-known  circuit  theory  definition.  Thus,  an 
admittance  or  impedance  descriptor  relating  v/(m)and  il(m)  defines  a  multiport  with  voltages 

v/(m)and  currents  il(m)  at  the  ports.  There  is  one  condition  that  has  to  be  satisfied  to  obtain  Fourier 

series  expansion  (4)  with  the  simple  multiport  power  definition  (6).  This  is  a  possibility  to  build 
a  set  of  orthogonal  basis  functions  (3)  on  the  surface  of  the  element.  How  to  proceed  with  such 
construction  is  not  obvious  for  a  polytope  in  a  general  case.  It  can  be  done  rigorously  only  for 
rectangular  or  brick  elements  that  have  rectangular  faces  and  adjoined  faces  belong  to  the 
orthogonal  planes.  Thus,  the  rest  of  this  paper  is  devoted  to  the  brick  elements  only. 

Basis  functions  for  a  rectangular  face 

Minimal  number  of  the  basis  functions  per  rectangular  face  for  a  3-diminsional  problem  is  two, 
to  account  two  polarizations  of  electric  field.  Each  basis  function  contains  one  constant  vector 
component  of  electric  field  and  one  constant  vector  component  of  magnetic  field.  The  basis 
functions  for  a  face  /  are  shown  in  Fig.  1.  Though  the  choice  of  the  directions  of  the  bases  is 
obvious,  it  is  not  so  with  the  magnitudes.  We  can  randomly  select  magnitude  of  the  electric  field 
and  define  magnetic  field  component  accounting  nonnalization  condition  (3).  In  general,  the 
basis  functions  may  be  expresses  as  follows: 
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where 

define 


e/(D  = 


=  «!*!>  /2/(l)=7L(«/Xe/(l))  =  A 


=  a2x2,  hl(2)  =  ^(n,  X el(2))  =  -p. 


(V) 


-1(2) 


Z°  =  \^/£  is  characteristic  impedance  of  the  medium.  There  are  multiple  possibilities  to 
a,.,/?,.,  /  =  1,2  ,  but  two  choices  correspond  to  some  physical  situations.  The  first  choice  is 


a, 

a2 


JL 

A2 ' 

,_7_ 

Al’ 


A  = 
A  = 


_7_ 

Al 

V_ 

A2 


where  coefficient  rj  has  to  be  chosen  to  satisfy  the  unit  power  nonnalization  condition  (3): 


(8) 

(9) 


The  ratio  of  magnitudes  of  the  electric  and  magnetic  fields  is  independent  of  r/ : 


7 


el(  1) 

.-*z 

7  — 

ei(2) 

\\) 

A2  °’ 

'(2) 

\(  2) 

-  —Z 
Al  ° 


(10) 


Z/(1)  and  Z/(2)  are  the  characteristic  impedances  of  the  parallel-plate  waveguides  with  electric 

walls  shown  by  thick  black  lines  and  with  magnetic  walls  shown  by  thin  dashed  lines  in  Fig.  1. 
With  the  choice  (8),  vl(m)  and  il(m)  define  normalized  voltage  and  current  of  the  parallel-plate 

waveguides.  With  77  =  1  the  admittance  descriptor  is  going  to  be  regular  admittance  matrix. 
Another  choice  of  the  coefficients  in  the  basis  functions  (7)  is  the  following: 


<7  /?[ 

^2  ~  Pi  ~ 


n 

VA1A2 

V 

VALA2 


(11) 


The  unit  power  transfer  nonnalization  condition  (3)  requires  the  same  value  (9)  for  the 
coefficient  77 .  The  ration  of  the  magnitudes  of  electric  and  magnetic  components  of  the  basis 
functions  with  (1 1)  is 


Z 


m 


eKD 

—  7  7  — 

eH  2) 

\  1) 

^0’  ^/(  2) 

h 

nl{  2) 

(12) 


This  definition  corresponds  to  the  basis  functions  defined  as  the  TEM  waves.  A  nonnalized 
admittance  matrix  obtained  with  either  definitions  of  the  surface  basis  functions  relates  surface 
expansion  coefficients  of  electric  and  magnetic  field  of  the  basis  functions  as  follows: 


/  =  V 


i.v  e  CN ,  T  e  C 


NxN 


(13) 


where  v  and  i  are  vectors  composed  of  unknown  surface  expansion  coefficients  (4),  N  is  the 
total  number  of  the  element  surface  basis  functions.  The  normalized  admittance  matrix  can  be 
formally  transformed  into  a  scattering  matrix.  Incident  and  reflected  waves  can  be  introduced 
formally  as 
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c+=\{v  +  i ),  c  =j(v-i ),  c+,c  eCN  (14) 

Corresponding  scattering  matrix  can  be  defined  as 

c-=S-c+,  SeCNxN  (15) 

It  can  be  expressed  through  the  admittance  matrix  as 

S  =  (U-Ynr)-(U  +  Yny  (16) 

where  U  is  the  unit  matrix  N  by  N  .  The  scattering  matrix  obtained  in  such  way  is  a 
generalized  scattering  matrix  and  does  not  depend  on  nonnalization  of  the  bases. 


Brick  element  descriptors 

Now  it  is  easy  to  show,  that  with  both  definitions  of  the  basis  functions  (7)  with  either 
coefficients  (8)  and  (9)  or  (1 1)  and  (9),  we  can  obtain  exactly  the  same  nonnalized  admittance 
matrix  of  a  brick  element  for  example.  The  admittance  matrix  of  a  brick  may  be  found  either 
with  formulas  for  its  elements  published  in  [2]  or  with  more  stable  generalized  matrix  procedure 
described  in  [4],  [5].  The  electric  field  components  in  [2]  have  been  chosen  as  constant  unit 
vector  functions  defined  on  the  faces  of  the  brick.  Magnetic  field  components  of  the  basis 
functions  in  [2]  have  been  also  chosen  to  be  constant  unit  vector  functions  defined  through  the 
electric  field  components.  With  the  projectors  defined  in  [2],  the  resulting  admittance  matrix  has 
to  be  normalized  in  general  case  to  be  converted  into  scattering  matrix.  Such  normalization  is  not 
required  in  case  of  the  cubic  element  that  is  investigated  in  the  paper  [2],  For  a  brick  element  the 
choice  of  basis  functions  and  projectors  of  paper  [2]  require  additional  normalization  of  the 
admittance  matrix  before  the  conversion  to  the  scattering  descriptor.  The  admittance  matrix  7 
defined  by  expression  (24)  of  paper  [2]  can  be  transformed  into  the  normalized  admittance 
matrix  (13)  as  follows: 

Ynr  —  Zo0  •  Y  ■  0  1  (17) 


where  0  is  a  diagonal  matrix  defined  as  follows: 

^AzAy  0  0  0  0  0 
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0  0  0  tJAxA v  0 

0  0  0  0  -JAxAv 


(18) 


Matrix  Ym  obtained  with  (17)  is  a  generalized  admittance  matrix  descriptor  of  the  element  and 
does  not  depend  on  the  choice  of  the  basis  functions.  It  is  interesting  to  investigate  the  un- 
nonnalized  admittance  matrix  with  the  two  different  sets  of  basis  functions  to  be  able  to  correlate 
it  with  the  actual  values  of  electric  and  magnetic  fields.  With  the  parallel-plate  waveguide 
definition  (8),  the  un-normalized  matrix  is 

Y  =  0,  -7 -0  1  (19) 

ppw  he  7/ 

where 
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Ax, 

Admittance  matrix  Y  relates  voltage  and  current  according  to  their  regular  definitions.  Wave 
channels  of  multiport  corresponding  to  brick  have  characteristic  impedances  defined  as 

Zc=Zo0,.0i*  (21) 

Normalization  of  the  matrix  Y  to  the  characteristic  impedances  of  the  channels  (21) 

Ynr  =  \[ZC  ■  Y  ■  ^Zc  gives  the  generalized  admittance  matrix  Ynr  defined  by  (17).  Admittance 
matrix  Y  can  be  also  normalized  to  characteristic  impedance  of  medium  Z0  or  to  50  Ohm 

before  transformation  to  the  scattering  matrix  (16).  It  does  not  mater  to  what  value  it  is 
normalized  as  long  as  it  is  appropriately  accounted  at  the  matrix  re-composition  stage  or  at  the 
interface  between  two  matrices.  With  the  TEM  waves  definition  (11),  the  un-normalized 
admittance  matrix  is 


Yfem  =  0  •  k  •  ©  1  (22) 

This  matrix  relates  directly  coefficients  of  expansion  of  electric  end  magnetic  fields. 
Characteristic  impedances  of  the  wave  channels  in  that  case  coincide  with  the  characteristic 
impedance  of  the  medium  Z0 ,  thus  normalization  to  Z0  produces  the  generalized  admittance 

matrix  Ynr .  Note,  that  all  descriptors  are  technically  equivalent.  The  normalization  does  not 

change  the  physics  of  the  problem  and  is  important  only  in  context  of  the  problem  at  the  element 
assembly  stage  or  to  construct  an  appropriate  energy  transition  procedure  between  the  elements 
in  time  domain  for  example. 


Conservativeness  of  the  brick  element 

To  illustrate  the  conservative  properties  of  the  brick  Trefftz  finite  elements,  let  us  consider  a 
cubic  element  with  Ax  =  Ay  =  Az  =  0.02  m  filled  with  vacuum.  The  element  size  is  approximately 
equal  to  wavelength  at  15  GHz.  Formulas  for  elements  of  admittance  matrix  from  [2]  have  been 
used  to  compute  element  descriptor.  The  admittance  descriptor  is  converted  into  generalized 
scattering  matrix  according  to  procedure  described  in  this  paper.  For  a  cubic  element  built  with 
the  collocation  projector,  the  scattering  matrix  has  only  one  distinct  element  with  non-zero 
constant  magnitude  0.5  and  linear  phase  -0.5kA .  Here  k  is  the  wave  propagation  constant  and 
A  is  the  element  size.  For  a  cubic  element  built  with  Galerkin  projectors,  there  are  just  two 
distinct  non-zero  elements  in  the  scattering  matrix  descriptor  [2].  Frequency  dependencies  of 
their  magnitudes  and  phases  are  plotted  in  Fig.  2  and  Fig.  3.  Note,  that  magnitudes  are  frequency 
dependent  in  that  case  and  phases  are  slightly  non-linear.  Both  cubic  elements  have  unitary 
scattering  matrix  descriptor  at  all  frequencies  that  can  be  proved  analytically.  Now,  let’s  build 
the  generalized  scattering  descriptors  for  a  brick  element  with  Ax  =  0.01  m.  Ay  =  Az  =  0.02  m  and 
filled  with  vacuum.  Frequency  dependencies  of  their  magnitudes  and  phases  for  a  brick  element 
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built  with  the  Galerkin  projectors  are  plotted  in  Fig.  4  and  Fig.  5.  The  magnitudes  values  are 
asymptotically  approaching  to  the  static  values  that  can  be  calculated  from  connection  of 
parallel-plate  waveguides.  If  scattering  matrix  is  not  generalized  the  asymptotic  is  going  to  be 
different  and  consistent  with  the  normalization  of  the  wave  channels.  The  scattering  matrix  is 
unitary  in  that  case  as  in  the  case  of  cubic  element.  Note  that  scattering  matrix  descriptor  of  the 
brick  element  built  with  the  collocation  projector  holds  this  property  only  approximately  if 
element  size  is  substantially  smaller  than  wavelength.  This  is  because  of  the  element  construction 
procedure  does  not  guarantee  the  conservativeness  with  collocation  projectors.  On  the  other  hand 
the  Galerkin  projectors  with  the  orthogonal  set  of  the  element  surface  basis  functions  allows  to 
enforce  the  conservativeness  of  the  element  descriptor. 
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Fig  1 .  Two  basis  functions  for  a  rectangular  face.  Black  solid  arrows  are  electric  field 
components  of  the  basis  functions.  Blue  dashed  arrows  are  magnetic  field  components  of  the 
basis  functions.  Red  arrows  are  vectors  normal  to  the  face. 
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Fig.  2.  Frequency  dependencies  of  magnitudes  of  non-zero  elements  of  scattering  matrix  of  the 
cubic  Trefftz  finite  element  of  Nikolskii. 


Fig.  3.  Frequency  dependencies  of  phases  of  non-zero  elements  of  scattering  matrix  of  the  cubic 
Trefftz  finite  element  of  Nikolskii. 
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Fig.  4.  Frequency  dependencies  of  magnitudes  of  non-zero  elements  of  scattering  matrix  of  the 
brick  Trefftz  finite  element  of  Nikolskii. 


Fig.  5.  Frequency  dependencies  of  phases  of  non-zero  elements  of  scattering  matrix  of  the  brick 
Trefftz  finite  element  of  Nikolskii. 
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